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ABSTRACT 

Recent monitoring programs of ultra cool field M and L dwarfs (low mass stars or 
brown dwarfs) have uncovered low amplitude photometric I-band variations which 
may be associated with an inhomogeneous distribution of photospheric condensates. 
Further evidence hints that this distribution may evolve on very short timescales, 
specifically of order a rotation period or less. In an attempt to study this behaviour in 
more detail, we have carried out a pilot program to monitor three L dwarfs in the near 
infrared where these objects are significantly brighter than at shorter wavelengths. We 
present a robust data analysis method for improving the precision and reliability of 
infrared photometry. No significant variability was detected in either the J or K' bands 
in 2M1439 and SDSS1203 above a peak-to-peak amplitude of 0.04 mag (0.08 mag for 
2M1112). The main limiting factor in achieving lower detection limits is suspected to 
be second order extinction effects in the Earth's atmosphere, on account of the very 
different colours of the target and reference stars. Suggestions are given for overcoming 
such effects which should improve the sensitivity and reliability of infrared variability 
searches. 

Key words: stars: low-mass, brown dwarfs - stars: variables: other - methods: data 
analysis - techniques: photometric - stars: individual: 2MASSW J1439284+192915 
- stars: individual: SDSSp J120358. 19+001550.33 - stars: individual: 2MASSW 
J1112257+354813 



1 INTRODUCTION 



Thermodynamical considerations show that dust particles 
form in the atmos phere s of ultra cool M, L and T dwarfs 
( e.g. A llard et al (kOOlT ); Burrows & Sharp ( |l999| ); Helling 
([200 1|); Lodders (|199S| ); Lunin e et al. ( |l989| ); Sharp & Hueb- 
ner (199C); Tsuji et al. (1996)). This has been confirmed in 



many cases from the modelling of their optical and near 
i nfrar ed spectral ene rgy d istributions (e.g. Legge t t et al. 
( ^00 1| ); Mar ley et al. (gOOf); Schweitzer et al. ( |200l[ )( |2002| )). 
Physical models suggest that as the effective temperature 
drops below about 2500 K, the photosphere becomes increas- 
ingly dusty due to increased gas condensation, whereas at 
even lower temperatures the dust opacity is reduced due 
to gravitational settlin g of t he dust below the photosphere 
(Ackerman & Marley ( ]200l| ); Tsuji ( |200l[ )). 

However, such models assume uniform horizontal distri- 
butions of this dust, whereas recent observational work has 
shown a significant fraction of ultra cool dwarfs to be pho- 
tometrically vari able a t the level of 0.01-0. 08 mag (Bailer- 
Jones & Mundt ([l999D(|200l[); Clarke et al. (]2002aj) (ppQ2b|) ; 



Email: calj@mpia-hd.mpg.de 



Gelino et al. ( |2002|) ; Marti n et a l. ( |200l| ); Nakajima et al. 
(bOQOl); Tinney & Tolley ( |l99S| )). Bailer- Jones & Mundt 
(2001) found evidence for a rapid evolution of features on 
the surfaces of a few field L dwarfs. Being field objects, they 
are presumably part of an older population which have long 
lost their disks (age > lOMyr), thus ruling out accretion hot 
spots or other star-disk activity as a cause of the variabil- 
ity. This leaves magnetic star spots or dust clouds as the 
most plausible candidates, although as yet there is no direct 
evidence for the physical nature of these surfac e features. 
Theoretical ar gume nts from Gelino et al. (2002) and Sub- 
hanjoy et al. ( [2002; ) argue against the existence of spots in 
ultra cool dwarfs due to the low ionization fractions. On the 
other hand, a non-monotonic variation with spectral type 
in the strength of the FeH band at 9896 A (which would 
be eff ected by iron condensation) noted by Burgasser et al. 

W ) 

may be indicative of an inhomogeneous dust cloud 
coverage, at least around the L/T transition. 

If the variability is due to dust, then understanding its 
physical, chemical and dynamical nature is central to deter- 
mining the fundamental physical parameters of ultra cool 
dwarfs (T c ff, abundances, ages etc.) and hence their for- 
mation mechanisms. Using the dusty and condensed dust 
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models of Allard et al. (2001), one of us has recently pre- 



dicted the observational signature s of d ifferent dust cloud 
and spot scenarios (Bailer- Jones (§00§)). This shows that 
there would be a distinct signature in the J and K bands for 
these phenomena, depending on effective temperature and 
feature size. 

In this paper, we present the results of a pilot project 
to look for variability in ultra cool dwarfs through near- 
simultaneous J and K' band monitoring. A comparison of 
the amplitudes and correlation of any variability can be com- 
pared with surface feature scenarios from various models. 



No correction was made for nonlinearity in the flux re- 
sponse of the pixels. In principle, a nonlinearity can intro- 
duce apparent variability if the stars do not remain exactly 
in the same place on the detector or if the stars or back- 
ground undergo a uniform change in brightness (e.g. due to 
changes in sky brightness or atmospheric extinction). How- 
ever, we have calculated that such effects introduce errors 
in the relative photometry of much less than 1% under good 
conditions, and no more than 1% in the very worst case. 
Note that bright reference stars were avoided to reduce the 
influence of nonlinearity. 



2 DATA ACQUISITION AND REDUCTION 

Details of the three targets monitored are given in Table [ij 
They were chosen from a list of sufficiently bright objects 
with appropriate RA/Dec for the observatory /season to rep- 
resent a range of L spectral types. The targets were ob- 
served over 16 nights in March/ April 2001 from the 1.23m 
telescope at Calar Alto, Spain, with the MAGIC camera, 
equipped with a 256x256 NICMOS3 array. The pixel scale 
was 1.15"/pix providing a field-of-view of 5'x5'. The ob- 
jects were observed in both the J (50% transmission points 
at 1.12 and 1.40/im) and K' (50% transmission points at 
1.93 and 2.27^im) filters. Each object was monitored contin- 
uously for several hours each night (alternating between J 
and K'), weather conditions permitting. Table @ shows the 
log of observations. The seeing varied between about 1.2" 
and 2.2" in J (poorer seeing images are not retained in the 
final light curves). Weather conditions varied across the run. 

A problem with the NICMOS arrays is the variation 
of the quantum efficiency (QE) across a single pixel (intra- 
pixel QE variation). On account of the coarse pixel scale and 
undersampled photometry, this could result in a variation in 
the detected flux of a non-variable star as its image moves 
relative to the center of a pixel. For this reason, each epoch 
(point in the final light curves) is derived from 45 frames 
taken with a dithering macro with non-integer pixel offsets: 
A field is observed at each of the nine points in a 3 x 3 square 
grid, with x/y distance 12" (10.4 pixels) between grid points. 
Five consecutive 5s integrations (repeats) are taken at each 
grid point, read out in double correlated readout (DCR) 
mode. The macro required 4.7 minutes to execute. Due to a 
variable bias level, the first repeat at each grid point could 
not be used, leaving a total of 36 usable frames. A clipped 
average of these points provides a single point in the final 
light curve (see section ^j). This procedure is necessary to 
achieve high precision photometry with infrared arrays. 

It was further found that the bias levels for repeats 2- 
5 were all different, although stable in time. Thus the flat 
fielding and sky subtraction had to be performed separately 
for each repeat. 

Flat fields were created on each night from the differ- 
ences between illuminated and non-illuminated dome flats. 
For each epoch, a sky image was then created from a me- 
dian combination of the nine frame positions for that repeat 
stage (each offset to have a common flux zero point to ac- 
commodate changes in the overall background) . This process 
removed all sources, leaving an image of the 2D background 
(predominantly sky brightness) variations which was sub- 
tracted from each frame. 



3 PHOTOMETRY 

A reference image was established and the positions of the 
target and reference stars in this propagated to all other im- 
ages of the field. This avoids shifting and interpolating im- 
ages, which can only degrade the data quality. Photometry 
was carried out using the CCDCAP aperture photometry 
program written by K. Mighell. CCDCAP was developed to 
do accurate aperture photometry with small, undersampled 
apertures, using a bilinear pixel inter polation algorithm (see 
appendix B of Mighell & Rich (1995)). After detailed inves- 



tigation of the magnitude error (see below) as a function of 
aperture size, aperture radii of 2.4 and 2.0 pixels were used 
for J and K' respectively. A 'hardness' of 1.0 was used in 
CCDCAP, which corresponds to the maximum subdivision 
of a pixel. A 'drift' of 1.5 pixels in CCDCAP allowed for 
optimal recentering of the aperture. 

Variability monitoring requires a reliable assessment of 
the photometric errors. Error determinations based on read- 
out noise and Poisson statistics in the source and sky gener- 
ally underestimate the true error because they ignore various 
systematic errors. The grid observing method outlined above 
is ideally suited for determining an empirical error measure 
based on the 36 frames taken in rapid succession. However, 
due to numerous bad pixels - up to four of the nine grid 
positions could be contaminated - simply taking the mean 
and standard error of these measures frequently gave poor 
results, because the bad pixels could bias the photometry 
by several magnitudes. As the PSF is undersampled, inter- 
polation over bad pixels (and cosmic rays) was not deemed 
reliable (and this can only be done to make reduction con- 
venient: it does not add information). Instead, an iterative 
rejection scheme was developed to remove bad photometry. 
First the median and its standard deviation, a, were de- 
termined from all 36 measures. If a exceeded a threshold, 
the highest 8 and lowest 8 measures were rejected, typically 
corresponding to all measures at four grid positions. The 
median is then re-established. Iterative clipping of points 
more than n-a about the median was then done to remove 
outliers (each iteration starting with all 36 measures so that 
re- inclusion was possible). The initial threshold stage is re- 
quired because the median of 36 points containing just a few 
bad measures can result in an unrepresentative median and 
hence large amounts of data being clipped. A threshold of 
0.2 mag and n=5 were used, although the results are not 
sensitive to these values. The iterative clipping converged 
after 0-3 iterations. This procedure was found to be very 
robust, both excluding obvious bad pixels but resulting in 



Limits on the infrared photometric monitoring of brown dwarfs 

Table 1. Targets observed. The reference is the discovery paper and provides the spectr al typ e 



and IAU name. The photometry is from 2MASS from the archive compiled by Kirkpatrick ( 2003 ) 



name 


IAU name 


SpT 


J 


K s 


reference 


2M1439 

SDSS1203 

2M1112 


2MASSW J1439284+192915 
SDSSp J120358.19+001550.33 
2MASSW J1112257+354813 


LI 
L3 
L4.5 


12.76 
14.02 
14.57 


11.58 
12.48 
12.69 


Reid et al. ( 
Fan et al. ( 
Kirkpatrick et al. ( 


2000 
2000 
2000 


) 
) 
) 



Table 2. Log of number of observations retained in final light curve for each object. 
AJD is defined as JD-2450000. 



Night 


AJD 


2M1112 


SDSS1203 


2M1439 


weather 






.1 


K' 


.1 


K' 


.1 


K' 




2 


1990.5 


24 


27 










some cloud 


3 


1991.5 


8 


10 










cloudy 


1 


1992.5 


12 


13 






41 


41 


mostly clear 


5 


1993.5 


11 


12 






21 


21 


clear 


6 


1994.5 










32 


32 


some cloud, later clear 


7 


1995.5 


3 


3 






12 


12 


some cloud, later clear 


8 


1996.5 














rain (no observations) 


9 


1997.5 










6 


6 


high humidity and cloud 


10 


1998.5 






7 


7 


26 


26 


mostly clear 


11 


1999.5 














fog (no observations) 


12 


2000.5 














fog (no observations) 


13 


2001.5 






39 


39 


10 


10 


clear 


11 


2002.5 






12 


42 


10 


10 


clear 


15 


2003.5 


13 


44 






6 


6 


clear 


16 


2004.5 










15 


15 


clear, high humidity 


17 


2005.5 






16 


16 


23 


23 


mostly clear 


Total 




101 


109 


104 


104 


232 


232 





more than 32 (of 36) measures being retained in the majority 
of cases. 

The mean of the remaining measurements was taken 
to be the final magnitude for that epoch. As these are all 
independent measurements of the same thing|] the appro- 
priate error, 8m, in this measure is the standard error in 
this mean, i.e. oj^fN . We found that the theoretical errors 
produced by IRAF's apphot task (based on source Poisson 
noise and detector read-out noise) - using the same aper- 
ture sizes - were as much as ten times smaller, i.e. grossly 
underestimated. This demonstrates the importance of using 
empirical error measures in infrared photometry, at least 
with NICMOS3 arrays. 

Photometry was obtained in this manner for the target 
and reference stars in all frames. Differential photo metry was 
performed as described in Bailer- Jones & Mundt (2001) to 



produce a light curve for the target relative to the average 
flux of the reference stars. This removes first order (wave- 
length independent) variations in the sky extinction between 
epochs.^] Light curves were similarly constructed for each ref- 
erence star relative to all the other reference stars. 6, 8 and 



1 This assumes that the star does not vary over 4.7 minutes. If 
it does, then this error is an overestimate of the true error. 

2 This is true provided the extinction is constant across the 5' 
image field size averaged across the 70 cos <5 arcmins on the sky 
through which the field moves during the 4.7 min macro execution 
time. 



5 reference stars were retained for 2M1112, SDSS1203 and 
2M1439 respectively. 

In what follows, light curve shall always mean the rel- 
ative magnitude light curve, i.e. the magnitude of the star 
relative to its particular set of reference stars. By reference 
level we mean the time series of the reference magitude for 
the L dwarf target, i.e. the magnitude formed from the av- 
erage of the fluxe s of t he target's reference stars (see Bailer- 
Jones & Mundt (bOOlh section 3.3). 



4 LIGHT CURVE ANALYSIS 

The light curves for the target star and reference stars were 
plotted and examined visually f or feat ures. The x 2 test as 



used by Bailer- Jones & Mundt (2001) was also applied to 



look for variability across the whole data set and within 
individual nights. When interpreted as probabilities of vari- 
ability on the assumption of Gaussian errors, these \ 2 values 
often indicated significant variability (p « 0.01) in both 
the target and some of the reference stars. Assuming that 
not all of these really are intrinsically variable at this level, 
this implies one or more of the following: 1. the photomet- 
ric errors have been underestimated; 2. the errors are non- 
Gaussian; 3. the relative photometry is not representative 
of the intrinsic brightness of the stars. Point (1) seems un- 
likely given the thorough testing of the method described in 
section |^ to evaluate the errors. Point (2) is generally true 
due to outliers: this will inflate all % 2 values for all stars, 
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meaning that one should adopt a conservative threshold for 
flagging variability (i.e. a small value of the probability, p, or 
equivalently a higher \ 2 value) . However, this would not be 
magnitude dependent so would not lead to a differential ef- 
fect between the target and reference stars. Thus the relative 
values of \ 2 can still be used as an indication of variability. 
Point (3) could be problematic on account of colour effects 
in the extinction variation, as will be discussed in section ^ 
As we have near-simultaneous observations in two 
bands, J and K', we introduce a parameter to look for cor- 
related changes in the relative magnitudes. For a given star 
at epoch number t, this parameter is defined as 



Q(t) = Y, 



mj(t') - mj(t' - 1) m K '(t') - m K ,(t' - 1) 



where mj(i') is the relative J magnitude of the star at 
epoch t' and similarly for m K i(t'). Smj is the error in 
mj(i') — mj(t' — 1), obtained from the quadrature sum of 
the empirical error measures, 8m, for each epoch (see sec- 
tion ^). Each of the two terms on the right hand side of this 
equation gives the change in relative magnitude from one 
epoch to the next in the units of the random error; signifi- 
cant changes have a modulus larger than unity. As the two 
terms, one for J and one for K', are multiplied together, cor- 
related changes in the two bands give a positive contribution 
to the sum; anticorrelated changes a negative contribution. 
Thus if the star shows a series of correlated changes, Q will 
become more positive; if they are anticorrelated it becomes 
more negative. If changes are random, i.e. sometimes cor- 
related and sometimes anticorrelated, Q will do a random 
walk; in particular it has an expectation value of zero and 
variance of t (the number of epochs). By comparing the vari- 
ation in Q for the target with that for its reference stars, we 
can see whether the target tends to show a greater level 
of colour correlated (or colour anticorrelated) changes than 
the reference stars. This is summarized by Q s , the difference 
between the maximum and minimum values of Q across all 
epochs. 

As cool spots and dust features effect the radiative flux 
over a wide wavelength range in the optical and infrared, 
the evolution of such features will produce some kind of 
correlated or anticorrelated change in the J and K' bands 
and hence some pattern in the Q parameter. This is the 
motivation for the Q parameter. The specific relative am- 
plitudes in these bands depends on the physical mechanism 
(spots or clouds or something else) and of course on the 
specific atmospheric models used. For example , mak ing pre- 
dictions based on the models by Allard et al. ( 2001 ), we see 
that dust variations in early L dwarfs produce anticorrelated 
J/K' variation s, wh ereas star spots produce correlated ones 
(Bailer- Jones (2002)). Different relative amplitudes may be 
predicted when making different assumptions in the models, 
e.g. with a different dust grain size distribution or treatment 
of convection. 

Our J and K' observations are not strictly simultane- 
ous. However, for SDSS1203 and 2M1439, a strict pairing of 
K' observations immediately after J observations was main- 
tained in the final light curves, and Q has been evaluated 
using this pairing. For various reasons this was not possible 
with 2M1112, so Q values have not been calculated for this 



Table 3. Photometry of targets and their reference stars. J r and 
K' r are the magnitudes of the reference stars and (J-K') r their 
colours, relative to the target (and averaged across all epochs). 
Sm d is the average across epochs of the relative mag nitude errors, 
as given by equation 4 in Bailer- Jones & Mundt j2001 ). Q s is 
defined in the text. 





Jr 


5m d 


K' r 


5m d 


(J-K')r 


Qs 


SDSS1203 


0.000 


0.016 


0.000 


0.012 


0.000 


9 


ref. 1 


-1.620 


0.010 


-0.849 


0.013 


-0.771 


7 


ref. 2 


0.578 


0.012 


0.448 


0.014 


-1.026 


15 


ref. 3 


0.829 


0.021 


1.757 


0.035 


-0.928 


11 


ref. 4 


0.640 


0.019 


1.614 


0.030 


-0.974 


7 


ref. 5 


0.932 


0.026 


1.701 


0.033 


-0.769 


9 


ref. 6 


1.467 


0.037 


2.562 


0.067 


-1.095 


7 


ref. 7 


1.242 


0.035 


1.516 


0.032 


-0.274 


6 


ref. 8 


1.433 


0.033 


2.376 


0.053 


-0.943 


10 


2M1439 


0.000 


0.016 


0.000 


0.018 


0.000 


15 


ref. 1 


-0.559 


0.017 


0.101 


0.023 


-0.660 


8 


ref. 2 


0.707 


0.018 


1.431 


0.024 


-0.724 


23 


ref. 3 


1.953 


0.026 


2.699 


0.042 


-0.746 


8 


ref. 4 


2.257 


0.034 


3.102 


0.058 


-0.845 


12 


ref. 5 


2.557 


0.038 


2.901 


0.050 


-0.344 


15 



field. A summary of the relative photometry, errors and Q s 
is given in Table ^. 



5 RESULTS 

SDSS1203 Calculated across the entire data set, SDSS1203 
gives x 2 values of 200 and 306 for J and K' respectively. 
These formally correspond to very small p values (a test 
with 103 degrees of freedom) but, as noted in section |i|, we 
may expect \ 2 to be large due to outliers across this large 
number of observations. Indeed, the reference stars also give 
large \ 2 values: ranging from 114 to 158 (for J) and 98 to 
155 (for K'). Ignoring the formal p values and just looking 
at the relative x 2 values, we see that SDSS1203 is indeed 
more variable than any of the reference stars in both filters. 

Fig. |l| shows the J and K' light curves across the four 
nights of observation. There are several points which deviate 
from the mean by more than 3a, although generally the 
variability is quite weak (< 0.05 mag deviation). On night 
13, there is a dip in K' at around AJD 2001.4. A similar 
but opposite effect is seen in reference star 2 (although not 
ref. 1) and the egress from this dip is accompanied by a rise 
in the reference level, so this is probably a telluric effect. 
The 4er jump in J at AJD 2001.51 is not accompanied by 
any similar feature in the reference stars or reference level so 
could be intrinsic to SDSS1203, although it is neither large 
nor that statistically significant. The J dip at AJD 2002.63, 
on the other hand, is accompanied by a dip in the reference 
level (but not the reference stars). This could be due to a 
colour dependent variation in the extinction (see section ^J). 
We also note that SDSS1203 appeared to be brighter in K' 
at the beginning of night 17, an effect which is marginally 
significant and cannot obviously be reduced to the behaviour 
of the reference stars or reference level. 

The Q parameter in Fig. [l] shows no long lasting colour 
correlated variations. Moreover, the range of Q, given by Q s , 
is no larger than for the reference stars (Table H). The J/K' 
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Figure 1. Light curves for SDSS1203, one plot for each night (nights 10, 13, 14 and 17). The open circles (blue) are 
for the J band and the closed triangles (red) for the K' band. The middle panel for each night shows the relative 
light curve of SDSS1203 (with a common magnitude zero point for all nights equal to the average relative magnitude 
across all epochs). The top panel for each night shows this light curve in units of the photometric errors at each 
epoch, which can be used to judge the significance of any deviations. The bottom panel plots the Q(t) value, defined 
by the equation in section |l| 



correlation coefficient, p, is slightly larger for SDSS1203 than 
the other stars (p=0.34 against p=0.07-0.24 for the reference 
stars), although it is still small. 

We conclude that there is no good evidence for variabil- 
ity in SDSS1203 from these data. SDSS1203 was fo und t o 
be variable in the I-band by Bailer- Jones & Mundt (2001), 
primarily due to an apparent brightening of the object over 
a duration of one to two hours. As far as we could tell, 
those data were taken under good observing conditions near 
culmination of the object (i.e. very small airmass changes) 
so we do not believe that result was an artifact due to the 
Earth's atmosphere. 

2M1439 As with SDSS1203, this object shows more signif- 
icant variability than any of the reference stars for both J 



and K' according to the \ 2 test applied to all data. However, 
its Q values, as well as the correlations between J and K' 
and J and J— K' are no larger than for some of the reference 
stars. 

It can be seen in Fig. | that the K' light curve of 2M1439 
dims by 0.14 mag over a period of five hours from AJD 
1992.4 on night 4. Although no such trend (or opposite 
trend) is seen in the light curve of the reference stars in- 
dividually, the opposite trend is seen in the reference level 
with a scale of 0.4 magnitudes, indicating that all of the ref- 
erence stars increased in brightness (Fig. ^|). This coincides 
with the field rising - from an airmass of 2.6 to 1.05 - so 
is presumably just due to the usual decrease in extinction 
with decreasing airmass. However, while the reference stars 
brighten by 0.4 mag, 2M1439 brightens by only 0.26 mag. 
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Figure 2. Light curves for 2M1439 (see caption to Fig. |l| for an explanation). Light curves are shown only for 6 of 
the 11 nights on which data were taken (see Table Eh, specifically for nights 4, 5, 6, 7, 10 and 12. 
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Figure 3. Relative magnitude light curves for 2M1439 (top row), the first (i.e. brightest) reference star (second row) and the second (i.e. 
next brightest) reference star (third row). The bottom row shows the variation in the reference level for 2M1439. The reference level is the 
magnitude formed from the average of the fluxes of all five reference stars shown in Table [}] (with the average of this level over all epochs 
subtracted). As in the other figures, the open circles (blue) are for the J band and the closed triangles (red) for the K' band. The left panel 
is for night 4, the right panel for night 7. 



This is possible if the extinction coefficient has a strong gra- 
dient across the K' band and if 2M1439 has a different flux 
gradient from the reference stars. In this case, second or- 
der (colour) effects in the terrestrial extinction may become 
important (see section |^) . The same effect is seen in the J 
band (reference level brightens by 0.35 mag, 2M1439 dims 
relative to this by 0.07 mag). 

However, on night 7 we see the same brightening of the 
reference level due to decreasing airmass (0.65 mag in both 
bands from airmass 2.8 to 1.05) but without any trend or sig- 
nificant variability in the light curves of either the reference 
stars or 2M1439 (see right hand panel in Fig. ^). Clearly, 
different atmospheric conditions must have prevailed on this 
night than on night 4. 

Correlated changes of up to 1 magnitude in J and K' 
are found in the reference level at the beginning of night 
6, caused by thin clouds and humidity. Correlated changes 
with amplitude up to 0.1 mag are also seen in the light 
curves of the reference stars and 2M1439, indicating that 
these large extinction variations have not cancelled out in 
the relative photometry. This is presumably caused by water 
vapour changes associated with the clouds (see section ^|). 

In conclusion, most of the features which account for 
the higher variability in 2M1439 (compared to the reference 
stars) appear to be related to terrestrial atmospheric phe- 



nomena. 



2M1439 has previou sly be en monitored in the I-b and b y 
Bailer- Jones & Mundt ( |200l| ) and by Gelino et al. ( fc002|) . 
Neither found good evidence for variability above about 
0.01 mag, the former from temporally-dense observations 
(48 observations spread over four nights) and the latter 
over longer timescales (a few observations per night over 40 
nights) . 

2M1112 This light curve shows some variations which at 



first sight would be of significance, but on closer inspection 
appear to be due to variations in the reference level. These 
are not qualitatively different from the "problematic" 
variations seen in the other two targets, discussed above. 

We performed a frequency-domain search for periodic vari- 
ations using two methods, the CLEAN and Lomb-Scargle 
periodograms. These methods yielded consistent results and 
showed no evidence for periodic variations in SDSS1203 or 
2M1439. As implemented, these methods are known to be ef- 
fective at dete cting even low signal-to-noise periodic signals 
(Lamm et al. (fc003|)). 



6 SECOND ORDER EXTINCTION EFFECTS 

The method of relative photometry employed in this pro- 
gram (and all those referred to in section nf) makes the im- 
plicit assumption that the Earth's atmospheric extinction 
is the same for the target star as for the reference stars. 
However, the extinction coefficient varies with wavelength, 
so if the target and reference stars have different spectral 
energy distributions (SEDs), their wavelength integrated ex- 
tinctions are also different. This is relevant when monitoring 
ultra cool dwarfs, as the reference stars will usually be sig- 
nificantly hotter and hence bluer in the optical and near 
infrared (see Table |^). This would not matter for rela- 
tive photometry if the extinction did not vary in time or 
space. But extinction does vary polychromatically with air- 
mass and atmospheric conditions, thus introducing a change 
in the relative magnitude between stars with different SEDs. 
Such effects are referred to as second order extinction. 

In the visible, the extinction is dominated by molecu- 
lar (Rayleigh) and aerosol scattering and generally shows a 
smooth decrease with increasing wavelength. In this case, 
colour effects are either small or a broad band colour term 
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can be used to correct relative photometry (also as a func- 
tion of airmass). In the near infrared, however, the extinc- 
tion is dominated by molecular absorption, particularly of 
H2O, CO2 and 03, which show much sharper variations with 
wavelength. Precipitable water vapour and ozone concentra- 
tions in particular can vary rapidly in time and space. Thus 
even under apparently good observing conditions, second or- 
der extinction effects can have a significant impact on high 
precision, differential, broad band infrared photometry. Note 
that terrestial water clouds and fog are dominated by large 
(>> 1 /im) water droplets, so these cause a scattering of vis- 
ible and infrared light with only a weak wavelength depen- 
dence (A~ Q , where < a < 1). Nonetheless, some absorp- 
tion also occurs in the liquid drops themselves which does 
show rapid changes in some wavelength ranges, although in 
the optical at least the colour dependence o f terre strial cloud 
extinc tion appears to be small (Honeycutt (1971); Serkowski 
(1970)). However, where there are clouds there is presum- 
ably also a high water vapour column density, with the result 
that clouds would contribute to second order extinction. 

Unfortuntately, we cannot quantify the effects of second 
order extinction on our data, because we have no indepen- 
dent measure of the atmospheric constituents - in particular 
the column density of the precipitable water vapour - and 
we would further require an appropriate atmospheric model 
matched to the Calar Alto conditions at the time. Moreover 
we do not know the SEDs of the reference stars: the broad 
band J— K' colour is too undersampled to determine the in- 
tegrated extinction from a water spectrum (Young (1989)). 

An idea of the scale of the problem can be found from 
reference to the models of J and K extinction at Kitt Peak 
for a cool giant (T cff =4000 K) and Vega (T cff =9650 K) from 
Manduca & Bell (1979). At low water column density, the 
differential extinction between the two stars (E4000— E9660) 
is 0.049 mag at J and 0.002 mag at K. Increasing the wa- 
ter column density by a factor of 35 changes these to 0.067 
mag at J and 0.001 mag at K. Thus increasing the precip- 
itable water density increases the differential extinction (so 
decreases the apparent brightness of a cool star relative to a 
hotter star) by approximately 0.02 mag at J and —0.001 mag 
at K0 The actual precipitable water column density expe- 
rienced during our observations could be larger, so it is not 
implausible that water vapour changes cause the apparent 
variability we see in the L dwarfs. Contrary to the figures 
given above, however, we generally see similar amplitude 
variations in J as in K'. But our K' filter extends further into 
the blueward telluric absorption band than does the John- 
son K band of Manduca & Bell and these figures are quite 
sensitive to the band profile. Also, our target (and probably 
reference) stars have much lower temperatures than 4000 K 
and 9650 K. 

Changes in airmass under stable conditions can also 
change the relative photometry. For example, Manduca & 
Bell show that the J band magnitude of a cool giant relative 
to Vega would be at least 0.01 mag higher when measured at 
airmass 2.0 than at airmass 1.0. The effect is much smaller 
in the K band, however, so this may not be cause of the 
dimming of 2M1439 on night 4 discussed in section |B| (see 



Fig. |^), as there the dimming is larger in K' than J. How- 
ever, this may again be due to the differences between the K 
and K' bands and the SEDs of the stars. Recall that no such 
dimming was seen on night 7 for the same airmass change, 
perhaps indicating different causes of extinction on these 
nights. 

The fact that both SDSS1203 and 2M1439 showed 
larger \ 2 values for J and K' than the reference stars can 
be attributed to the likelihood that the reference stars have 
more similar SEDs to each other than to L dwarfs, so that 
differential extinction effects between any one reference star 
and all the other reference stars (which form its reference 
level) are much smaller. Hence the reference stars show less 
variability. 

Second order extinction is less of a problem in the op- 
tical (A < 0.65 /im), but here most brown dwarfs are far 
too faint to be monitored. The I band (~ 0.78-0.92 /mi) is 
a good compromise. This too is intersected by a water ab- 
sorption band centered at 0.81 /im, but this band is much 
weaker t han th ose effecting the J and K bands (e.g. McCord 
& Clark (]1979|)). 



3 Derived from Tables Ha and lib of Manduca I 
niters at an airmass of 1.5. 



Bell for Johnson 



7 CONCLUSIONS 

The prime requirements of any search for variability are (1) 
a measurement of the scale of the variability intrinsic to the 
source, and (2) an accurate determination of the observa- 
tional errors in these measurements. However the analysis 
is done, a variability detection relies on the former being 
significantly larger than the latter. 

We have developed an observation and data reduction 
technique which gives an accurate and reliable determina- 
tion of the photometric errors, arising from the source, back- 
ground, detector and effect of the data processing (see sec- 
tion |). 

One of the main problems with ground-based monitor- 
ing is variations in the extinction (scattering plus absorp- 
tion) of the Earth's atmosphere. The standard procedure of 
differential photometry was used to give a measure of in- 
trinsic variability. Using this, some evidence for variability 
was found in the observed L dwarfs. However, closer analy- 
sis has shown that this measure is probably "contaminated" 
by second order extinction effects in the Earth's atmosphere 
(see section ^) . Although we have no direct evidence for this, 
it seems the most plausible cause of the observed variabil- 
ity given the features we have described and the elimina- 
tion of other potential causes. We therefore conclude that 
we have no good evidence for intrinsic variability in any of 
the three L dwarfs monitored. None of these objects showed 
significant correlations in their J and K' light curves, as ev- 
idenced by correlation plots and our Q parameter. Some 
such (anti)correlation would be expected if the variability is 
caused by cool spots or dust features. Upper limits on the 
peak-to-peak amplitude of persistent variability are set by 
the scatter in the light curves under the most stable Earth 
atmospheric conditions. For both J and K', these limits are 
set at 0.04 mag for SDSS1203 and 2M1439 and 0.08 for 
2M1112. 

As late M, L and T dwarfs are at their brightest in 
the near infrared (0.9-2.5 /im) and show spectral signatures 
which could discriminate between various surface feature 
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scenarios (Bailer- Jones (2002)), this is a desirable wave- 



length range for variability monitoring. We have shown that 
high precision differential photometry (random errors of less 
than two percent) can be achieved on L dwarfs with infrared 
arrays. However, as a measure of the intrinsic variability in 
L dwarfs, broad band differential infrared photometry ap- 
pears to be limited in accuracy to a few percent by variations 
in terrestrial molecular extinction, in particular precipitable 
water vapour. These second order extinction effects could be 
reduced in one or more ways: 

(i) use of passbands which avoid molecular extinction, in 
particular use of better designed (narrower) J and K' filters 
to avoid the strong H2O bands which generally intersect 
these "standard" filters^; 

(ii) use of narrower band filters, which reduces the effect 
of second order extinction^] (this would obviously require a 
larger telescope or longer integration times; the latter will 
be unacceptable for detecting short-term variability); 

(iii) use of specially designed filters to monitor the dom- 
inant molecular absorption (in particular H2O) from which 
the broad band relativ e pho tometry could, in principle, 
be corrected (Angione (1987)). This could be done using 
a multi-channel camera. However, additional spectroscopic 
observations are required to calibrate the measures and 
so apply the correction. Nonetheless, some kind of time- 
dependent monitoring of extinction would at least allow one 
to recognise when observing conditions were variable; 

(iv) observing from very dry sites under atmospherically 
stable conditions. Space-based observations remove all prob- 
lems related to second order extinction. A strong discrim- 
inant between candidate surface features in brown dwarfs 
is the variability signatures in the water absorption band s 
at 1.35-1.45/im and 1.80-2.00/wn (e.g. Bailer- Jones ( fc002|) . 
These bands cannot be reliably monitored from the ground, 
providing another argument for space-based observations. 
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